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Abstract 

An uplink system with a single antenna transmitter and a single receiver with a large number of 
antennas is considered. We propose an energy-detection-based single-shot noncoherent communication 
scheme which does not use the instantaneous channel state information (CSI), but rather only the 
knowledge of the channel statistics. The suggested system uses a transmitter that modulates information 
on the power of the symbols, and a receiver which measures only the average energy across the antennas. 

We propose constellation designs which are asymptotically optimal with respect to symbol error rate 
(SER) with an increasing number of antennas, for any finite signal to noise ratio (SNR) at the receiver, 
under different assumptions on the availability of CSI statistics (exact channel fading distribution or 
the first few moments of the channel fading distribution). We also consider the case of imperfect 
knowledge of the channel statistics and describe in detail the case when there is a bounded uncertainty 
on the moments of the fading distribution. We present numerical results on the SER performance 
achieved by these designs in typical scenarios and find that they may outperform existing noncoherent 
constellations, e.g., conventional Amplitude Shift Keying (ASK), and pilot-based schemes, e.g.. Pulse 
Amplitude Modulation (PAM). We also observe that an optimized constellation for a specific channel 
distribution makes it very sensitive to uncertainties in the channel statistics. In particular, constellation 
designs based on optimistic channel conditions could lead to significant performance degradation in 
terms of the achieved symbol error rates. 
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I. Introduction 

Accurate channel state aequisition is eritieal to realizing many of the diversity and beam¬ 
forming benefits in existing high speed eommunieation systems. With an inereasing number of 
antennas and operating bandwidths, the eomplexity and overhead of ehannel state aequisition 
grows proportionally. However, in a reeent result [1], [2] for a narrowband massive SIMO system 
with independent and identieally distributed ehannel eoeffieients at eaeh of the reeeiver antennas, 
it was shown that one ean aehieve a sealing law in the number of reeeive antennas whieh is 
no different from that with perfeet ehannel state information (CSI). The aehievability uses only 
an average energy measurement aeross all the reeeiver antennas. In another reeent work [3], it 
was shown that employing a similar energy-based seheme in a noneoherent massive wideband 
SIMO system leads to optimal eapaeity sealing when both the number of reeeive antennas and 
the bandwidth go to infinity jointly. These asymptotie eharacterizations suggest that when the 
number of antennas is large, the requirements on the ehannel state information aequisition may 
not be very stringent. 

Motivated by these asymptotie results, in this work, we eonsider a finite number of reeeive 
antennas in a massive SIMO system and address the problem of optimizing the transmit 
eonstellation points based on an asymptotieally tight and analytieally traetable upper bound 
on the Bit Error Rate (BER) for any given SNR, under different assumptions on the knowledge 
of the CSI statisties. We also present robust eonstellation designs given imperfeet knowledge of 
the ehannel statisties. Keeping in mind the faet that the ehannel moments ean be estimated on 
time seales larger than those needed for estimating instantaneous phase (for whieh we typieally 
need resouree-eonsuming training sequenees), in our work we investigate how many antennas 
we need for a eertain BER performanee. Through analysis and simulations, we find that the 
suggested sehemes ean outperform several existing noneoherent sehemes, e.g., eonventional ASK 
modulation, as well as pilot-based sehemes, such as PAM modulation when the eoherenee time 
of the ehannel is small. 

To the best of our knowledge, this line of work, with initial results appearing in [4], is the first 
to eonsider an energy-based eneoding and deeoding proeedure for a noneoherent large antenna 
system and demonstrate its performanee gains over other noneoherent and eoherent sehemes. 
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A. Prior Work 

Noncoherent eommunieation, i.e., eommunieation without instantaneous ehannel information, 
has attraeted a lot of researeh due to either the diffieulty of ehannel state aequisition, or the 
need for low hardware eomplexity and low energy eonsumption. The earliest ineamations of 
noneoherent systems were mostly motivated by the simplieity of the reeeiver eireuitry: the use 
of envelope deteetors ean be traeed baek to the well-studied quadrature and square law reeeivers 
[5], [6] employed in the noneoherent deteetion of several well-known modulation sehemes, 
sueh as Frequeney-Shift Keying (FSK), Amplitude Shift Keying (ASK) [7] and Pulse Position 
Modulation (PPM) [8]. However, the speetral effieieney of these systems is inferior to that of 
their eoherent eounterparts. Henee, with the paueity of speetrum resourees in eellular and with 
the advanees in deviee manufaeturing, eoherent systems using phase aequisition eireuitry at the 
reeeivers mostly replaeed noneoherent systems. 

However, as the demand for mobile data in wireless broadband eommunieations inereases 
dramatieally every year, there is an evident trend towards higher and higher earrier frequeneies 
and large antenna arrays [9], [10]. In these seenarios, the issues of simple eireuit designs, 
inexpensive hardware eomponents and energy effieieney become as crucial to system design 
as speetral effieieney [9], [11], [12]. While the number of RF ehains goes up with an inereasing 
number of antennas, thereby eausing inereased eomplexity and energy eonsumption, hardware 
impairments sueh as phase noise and I/Q imbalanees also beeome more severe at both the 
transmitter and the reeeiver. Arehiteetures whieh use simple, robust and energy effieient designs 
are thus important to realizing many of the performanee benefits of massive MIMO. 

Initial investigations on the rate loss ineurred by even a small training overhead in a massive 
SIMO system show that in several seenarios of low SNR, high mobility, or a large line of sight 
(LOS) eomponent, a noneoherent system aehieves a better probability of bit error than a eoherent 
system for the same effeetive rate [13]. Thus, noncoherent sehemes seem to be a eompetitive 
alternative to eoherent sehemes in the design of praetieal large-antenna systems. 

Within this line of researeh, [14] foeuses on the noneoherent ML deeoder and proposes signal 
eonstellation designs using a metrie motivated by a union bound on the probability of error in 
the high SNR regime. Similar metries, also motivated by a high SNR analysis, are presented in 
[15] where the worst-ease ehordal distanee is employed to plaee the eodewords as far apart as 
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possible. An alternative to the ML detection of the transmitted symbols is considered in [16], 
where the authors consider the problem of joint channel and transmitted symbol estimation, 
and propose a minimum distance criterion for code design using the Generalized Likelihood 
Ratio Test (GLRT) in the AWGN channel at high SNR. The GLRT is necessary when the 
channel statistics are not known precisely and it involves maximizing the noncoherent likelihood 
function over all possible values of the channel statistics. While these decoders are applicable 
to a very general class of channels, in this work, our focus is on understanding and optimizing 
the performance of energy based decoders with an increasing number of antennas. 

The rest of the paper is organized as follows. We present the system model in Section II, and 
summarize relevant work on its asymptotic characterization in Section III. Then, Section IV-A 
presents the constellation design problem and Section IV-B describes the solution to this problem 
when the channel distribution is perfectly known. Section IV-C shows how the constellation 
design problem can be done when only the first few moments of the fading distribution are 
known, and Section IV-D addresses the case when even the moments are not known perfectly. 
Finally, in Section V, we present plots showing the performance of the suggested schemes with 
representative statistics. Section VI summarizes this work. 

B. Notation 

We use [k] to denote the set {1,2, • • ■ , k} where k is an integer. is the set of all complex¬ 
valued matrices of size nx m. For a matrix H G the (i, j)-th element is denoted by Hi^ 

and for a vector h G the i-th element is denoted as hi. Re(-) and Im(-) represent the real 
and imaginary terms, respectively. CJ\f{n, R) represents the distribution of circularly symmetric 
complex Gaussian (CSCG) random vectors with mean vector fx and a covariance matrix R. The 
symbol = is used to denote a definition. V refers to a set of power levels that the transmitter 
uses and n is used to denote the number of receive antennas. We refer to a {K, 7) Rician fading 
channel as a channel with Rician fading (iF-factor in dB units and unit second moment) and 
additive Gaussian noise with power = —7 in dB. 
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II. System Model 

Consider one single antenna transmitter in a flat fading channel and a receiver with n antennas, 
where n is a large (but finite) number. The system is represented as 

y = hx + V, (1) 

with y e X e C, V G h e and each x* ~ CA/’(0,(T^), hi ~ f{h), such that 

E[/rj] = fj,, E[|/rj —/ip] = a^, and f{h) is the probability density function of the channel 
distribution. For normalization purposes and for notational simplicity, we also assume that 
E[|/rip] = 1 and E[|xp] = 1 so that parameters such as long-term shadowing, path-loss and 
antenna gain are incorporated in the Then, the average SNR per antenna at the receiver for 
this model is 

A E[|/r,p] 1 

h = -^ 

We further assume that the density function f{h) is such that, for any fixed x G C, the moment 
generating function of |//ip, i.e., exists and is twice differentiable in an interval around 

9 = 0. Many fading distributions fall within this model, e.g., Rayleigh and Rician fading [17], 
in which case hi ~ CN‘{^, a^). 

Note that an important aspect of this system model, similar to many works in the massive 
MIMO literature, is the assumption that the channel realizations across the antennas are i.i.d. 
random variables. While this assumption is not typically accurate in practice, some recent 
measurements [18], [19], [20] suggest that, despite the statistical difference between the measured 
multiple antenna channels and the i.i.d. channels, most of the theoretical conclusions made under 
the independence assumption are still valid in real massive MIMO channels. 

Motivated by our recent asymptotic results in [1], [2], [3], this work focuses on symbol-by- 
symbol encoding and decoding schemes that use an average energy-based transmitter and receiver 
design. This means that information is modulated in the power of the transmitted symbols, |xp, 
and the receiver estimates only the average power of the received signal, We describe this 
next. 


July 20, 2015 


DRAFT 



6 


A. Transmitter architecture: energy encoder 

The transmitter encodes information only in the power of the transmitted symbols, i.e., it 
transmits symbols with power levels from a codebook 

'P = {Pi,P2,-" ,Pl}, 

where pk G M+, subject to an average power constraint 

1 ^ 

< 1 , 

k=l 

assuming equiprobable signaling. Here pk ^ V is the power level of the symbol and L is 
the cardinality of V. Note here that we do not encode information in the phase. 


B. Receiver architecture: energy decoder 

Assume the user transmits a symbol whose power is the constellation point from V, i.e, 
Pk- In order for the receiver to detect pk, it only computes the following statistic 

W ^ ELi \y^ ^ (2) 

n n ’ 

i.e., it estimates only the average received power across all its antennas. Based on its knowledge 
of the statistics of the channel, the receiver divides the positive real line into non-intersecting 
intervals or decoding regions 

V = 


(each Xk corresponding to each pk G "P), and returns 

(3) 

The probability of error Pe{pk) when the power level G P is transmitted, and the 
average Symbol Error Rate (SER) Pg for any fixed constellation size L are defined as 

1 ^ 

Pe{Pk) = Pr{k^k}, P. --5^Pe(Pfc), (4) 

^ k=l 

respectively, assuming equiprobable signaling. 
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C. Discussion 


The use of energy-deteetion-based transmission and decoding is motivated by the fact, proved 
in [2], that such an encoding and decoding method achieves the same BER as a noncoherent 
maximum likelihood (ML) scheme in the Rayleigh fading channel. To see this, as explained in 
more details in [2], assume the transmitter sends a symbol x G C. The log likelihood function 
for a noncoherent Rician fading channel, i.e., hi ~ CAf{^,al), is 


iog/r(y) = 


||y-/ra;l|p 


and therefore, the noncoherent ML decoder is 


+ nlog {Ti{a + al\x\^)) , 


k = argmax^^l^p^p^ vfc log (5) 

For /i = 0, i.e., Rayleigh fading, the noncoherent ML decoder depends only on ||y|p. Thus 
for suitably chosen decoding regions V, it performs as well as the ML decoder. In general for 
/r 7 ^ 0, energy based detectors are not optimal. However, as shown in Figure 3 for representative 
values for /r, the SER performance gap is typically small. 

The proposed architecture requires a very simple one-dimensional statistic of the received 
signals, which allows for a simplified RF chain design. A general noncoherent ML or coherent 
detector, on the other hand, requires much more complicated circuits. 


III. Constellation Design Optimization 


We are interested in the problem of minimizing SER for any fixed constellation size L and 
fixed n, i.e.. 


1 


= argmin —\og{Ps) 
n 

1 ^ 

subject to — E Pk < 1, 


( 6 ) 


fc=i 


0 < Pk, V/c e [n]. 

The optimization above is over all codebooks V and decoding regions V. Unfortunately, this is 
in general a difficult problem to solve. The scope of this work is to solve a specific relaxation 
of this problem motivated by the large n asymptotics (n is large but finite) for any channel 
distribution that has a m.g.f. Specifically, we consider maximizing the error exponent of SER 
with respect to n, or a second-order approximation of it. Before describing the relaxation we 
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first define the notion of the error exponent of SER with respect to n. Using this asymptotic 
characterization not only helps us make the constellation design problem tractable, but also gives 
us designs which work for little prior statistical characterization and decouples the effect of the 
number of receive antennas from the effect of the channel distribution. 


A. Error exponents 

Fix any codebook V. Define the receiver’s constellation points r{pk) to be the value of the 
average received energy when the transmitter sends the power level, i.e., 

r{Pk) =Pk + cr^- 


We now observe that 


||y||2 ||h^ + v||2 ||h||2 ||v||2 ^Re(h*v) 

-=-^-= - Pk H-h 2- 

n n n n n 

so, in the limit of large n, due to the law of large numbers and the independence of h and v, 

it follows that 



n^oo Tt 


For finite n, the statistic would deviate from the value r{pk). To characterize this deviation, 
it is helpful to define 




\hiy/^ + Oi 


E 


\hi^/^ + Ui 


\hi^/^+Vi\^ -r{pk) 


(7) 


as the random variation of the received energy at the antenna around its expected value 
r{pk). Note that {uk,i}i=i are independent realizations of the same zero-mean random variable 
Uk ~ Qkiu) whose m.g.f. is 


Mk{e) = 


( 8 ) 


In the above, Mk{9) depends on the statistics of the channel and the noise, and the power level 

Pk- 

In [2], we derived an upper bound for the objective in (6) as follows: 

1 ^ 

Ps< - E ^Q-nlR,k(dR,k) _|_ Q-nlL,k{dL,k)^ 

^ k=l 
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where 


lL,k{d) = sup {9d - \og{Mk{-9))) , lR,k{d) = sup {9d - \og{Mk{9))) , (10) 

e>o e>o 

are defined as the left and right rate functions of pk respectively. In addition, sup^^g^/(x), is 
the supremum of f{x) in A. The decoding regions are chosen as 


= {r{pk) - dL,k, ripk) + dn^k] , 


for some di^k and d^^k- Define 

Ik = min {lL,k{dL,k), lR,k{dji^k)) 


to be the rate function of the constellation point pk. Then, it was shown in [2] that 

j A l0g(P.) . . .... 

If, — hm-= minii., (11) 

n-5>oo n k&[L] 

i.e., the error exponent of SER, denoted as Jg, is the same as the worst rate function of the 
constellation points. In other words, for n large enough, the probability of error performance is 
dominated by the constellation point with the worst rate function. Therefore, the constellation 
points V and the corresponding decoding regions V could be chosen in such a way as to 
maximize the error exponent of SER, i.e., (error exponent maximization problem) 


maximize Je 
1 ^ 

subject to — E Pfc < 1, 0 < Pk- 

^ k=l 


( 12 ) 


Observe that problem (12) is a relaxation of problem (6) for any finite n, since its objective 
function is an upper bound on the objective function of the latter. Yet, note that due to (11) this 
relaxation would provide an asymptotically optimal constellation with increasing n, even if it 
does not explicitly solve the SER minimization problem exactly for any finite n (6). Eurthermore, 
an interesting aspect of (12) is that it depends only on the channel distribution, and not on the 
number of antennas. By decoupling number of antennas from the channel distribution it allows 
us to study achievable SERs with large n without explicitly solving (6). 

In [2] we showed the following about the left and right rate functions for any pk'. 


Lemma 1. The left and right rate functions lR,k{d), iR^kid), respectively, of the power level pk 
enjoy the following properties: 
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They satisfy 

Ito = lim ’-iff - 

d —>-0 (it d —>-0 (X 

with h ~ f{h) and v ~ CJ\f{0, 


2E|C/|]’ 


where Uk = \hk/pk + 'Wp - Pk - er^y 


They are non-negative, convex and monotonically increasing for positive d for a fixed non¬ 
negative Pk, and monotonically decreasing for non-negative pk for a fixed positive d. 

It holds that Iipif) = -^R,fc(0) = 0/or any non-negative pk. 


The above lemma provides important insights into the dependence of the rate functions on 
the power levels and channel statistics. Specifically, for a small d, which practically means 
large constellations, we can approximate the rate function using only the first few moments of 
the channel distribution (E[Ul] depends on the first, second, and forth moment of the fading 
distribution as we explicitly show in Section IV-C). Also, for a fixed d, increasing pk leads to 
smaller rate functions, i.e., worse SER performance, or in other words, the constellation points 
that correspond to high power levels generally experience worse performance than those with 
low power levels for decoding regions of the same size. We now describe how these analytical 
properties could be used to solve for the maximum error exponent as defined in problem (12). 


IV. Constellation designs 

A. Overview 

In this section we consider three cases of constellation designs to maximize the error exponent 
in (12), where each case corresponds to a different assumption on the availability of statistical 
information about the channel. 

• Case 1: Subsection IV-B presents a design which assumes that the encoder and decoder 
know perfectly the channel distribution. This constellation is denoted as Vk,-^- 

• Case 2: Subsection IV-C presents a design in which only the first four moments of the 
channel distribution, are perfectly known. We will show later that these four moments 
suffice for a near-optimal design at high rates. This constellation is denoted as Vk,'^ 

• Case 3: Subsection IV-D presents a design in which even the latter are imperfectly known. 
The corresponding constellation is where a is the uncertainty in dB around the nominal 
values K and SNR. 
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Algorithm 1: Bisection algorithm 
function: [t* ,V*] = Bisection( ) 

ti 0 , tji oo 

repeat 

f _ ii -\-tu 

^ ~ 2 

[V, St] = ConstellationDesign( t ) 

If {St < 1): ti = t; else = t 
until \tu — ti\ < e and \St — 1| < e 
return V* =V;t*= t; 


In addition to the constellation designs above, we also consider a minimum distance constellation 
design, denoted as that was proposed in [2], which has 

2(^-1) ^ , rri 

Pk = ^ ^ , for k e [L] 


with decoding regions 


Xn = 


0 , 


L-l 


2k-3 2k-1 
L- 1 ’ L-l 


Wke{2,--- Il = 


2 - 


L-l 



We now describe each of the above constructions. 


B. Perfect knowledge of channel distribution 

We first discuss the constellation design with perfect knowledge of the channel distribution at 
the receiver. Since the exact channel distribution is known, Mk{6) is also known at the receiver 
and transmitter for any chosen p^. Then, (12) can be written as 

maximize min {hAd^k), lR,k{dR,k)) 

{Pki(^L,k^dR,k}kG[L] kG[L] 

0 <Pl < P2 < ■ ■ ■ < Pl, 

(13) 

dL,k > 0 , dji^k > 0 , 

k=l 

assuming decoding regions of the form = {r{pk) — di^k ^^{Pk) + dR^k]-, where for simplicity 
we assume that dip = dR^i = oo. 
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Algorithm 2 describes in detail how to get the solution of the optimization problem (13) and a 
detailed proof is presented in Appendix A. To exemplify the procedure and provide an intuitive 
argument for the validity of the suggested construction we consider the case with L = 4 shown 
in Figure 1. The design is based on the following two properties that follow from Lemma 1: 

1) Both lL,k{d) and Ir, k{d) are non-negative and monotonically increasing functions of d for 
a fixed pk- This means that increasing the size of the decoding regions always helps to 
increase the resulting rate functions, and therefore increase the minimum amongst them. 

2) Both lL,k{d) and dn, k{d) are monotonically decreasing functions of pk for a fixed d. 

^R,l ^L,2 ^R,2 ^L,3 ^R,3 ^L,4 

■ ■ C2 ■ . 

Pi+o^ P2+o^ Pa+O^ P4+02 

Fig. 1; Example of the decoding regions using L = 4 

Based on these two properties, we have the following sequential construction: assume there 
exists a constellation with error exponent t* that satisfies the power constraint. This means that 
the left and right rate functions of all the constellation points at the receiver are at least t*. To 
find this constellation choose first the minimum possible value for pi. Then, choose the boundary 
of the decoding region to the right of r(pi) = pi + i.e., ci, as show in Figure 1, such that 

the right rate function of r(pi), /r,i(ci — r(pi)), is at least t* on the boundary. Then, choose 
the smallest p 2 such that r{p 2 ) > ci and the left right rate function of r{p 2 ), /l, 2 (^"(P 2 ) — ci), 
is at least t*. Note that choosing a higher p 2 is always an option but this will lead to a design 
that uses more power than necessary. We perform this procedure sequentially until we find pl. 
Then we check if the average power constraint is satisfied. If that is the case, the assumption 
that there exists a constellation with error exponent at least t* that satisfies the power constraint 
was correct. If not, we should discard this constellation, decrease t* and repeat the procedure. 
Proof that this procedure gives the optimal error exponent is presented in Appendix A. 

For L = 4, Figures 2-(a) and 2-(b) show the normalized empirical histogram of the received 
statistic in a Rayleigh fading channel of the suggested constellation: Vk,'^ and The 

circles and diamonds on the x-axis show the r{pk) and the Ck (boundaries of the decoding 
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n 



Fig. 2: Example of the histogram of for n = 100 antennas for Vk,^ and 7?™*" 


regions) respectively. Observe that for the constellation there is a significant overlap in the 
histogram and therefore the receiver experiences significant symbol error rates. Also, observe 
that as pk increases, the variation around r{pk) also increases, due to the special nature of the 
energy detector at the receiver. 

Furthermore, Figure 3 shows a numerical comparison of the SER of two systems that transmit 
using the same codebook V in Rician fading, but with different decoders; the first system 
uses the decoder in (3) and the decoding regions resulting from the above procedure, and the 
second system uses the ML noncoherent decoder (5). Observe that there is no difference in 
the performance of the two systems in Rayleigh fading. Moreover, even in Rician fading with 
a relatively strong LOS component, the difference in performance is small. This suggests that 
energy decoders can roughly match the performance of optimal decoder structures even when 
the fading distribution has a non-trivial LOS component. 

C. Perfect knowledge up to the forth channel moment 

The constellation design presented above assumes that the receiver knows exactly the channel 
statistics. This may not be realistic in a practical scenario. In this section, we relax this assumption 
and consider the scenario in which the encoder and decoders only know the first few moments of 
the channel distribution, i.e., up to the fourth channel moment. The latter is especially important 
since determining the exact small-scale fading channel models in some cases, such as millimeter 
wave frequencies, is still ongoing research, and may not be reliably known beyond the first 
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(a) Comparison of Energy and ML decoder in Rician fading with 
7 = 10 dB 



(b) /S' = 0 dB, 7 = 5 dB 

Fig. 3: (a) Comparison of the energy decoder with the noncoherent ML decoder in Rician fading, (b) Comparison 
of 'Pk,j as a function of the constellation size L in Rician fading 


few moments. Moreover, the first few moments of the underlying ehannel distribution eould 
potentially be estimated more effieiently (i.e., on large time seales) without the need for resouree- 
eonsuming training sequenees. 

To see why the knowledge of the first four moments is enough to design good eonstellations, 
let h = hre + jhim, where h^e, him ^ Using Lemma 1 leads to 


lR,k{dR,k) ~ lR,k{dR,k) — 


dl 


R,k 


2s{pk 


lL,k{,dk^k') ~ lL,k{dij^k') — 


dr 


L,k 


2s{pk 


(14) 


for small dR^k and di^k, with s{pk) = E[U|] = aipl + a 2 Pk + as, and 


«1 = ^[ Ke ] + E[/iL] + 2 - 1, 02 = 03 = 


These expressions follow from the Gaussianity of the noise and the faet that the noise and 
the ehannel are independent and that the noise is zero mean. Observe that this approximation 
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Algorithm 2: Constellation design: Perfect channel distribution knowledge 
[t*^V*] = Bisection( ); 


function; V = ConstellationDesign( t) 

Pi = 0; dL.i = oo; = oo; 
for fc = 1 , 2,3, • • • , L do 

if /c = 1 then 

Pi = 0 

else 

Find the smallest pk > Pk-i + cJit,/c-i such that lL,k{Pk - Pk-i - dn^k-i) = t 
dL,k = Pk — Pk -1 — dfi^k-l 

end if 

if k ^ L then 

Find the smallest d^^k > 0 such that lR^k{dR,k) = t 

end if 

= [pk + 0-2 - dL,k,Pk + 0-2 + dR^k] 

end for 

P = bl,P2,--- ,PL\,St = i^t=lPk ’ = {d^k}k=l 

return V,T>,St 

end ConstellationDesign 


depends only on the first, second and fourth moment of the channel distribution. For example, 
in the case of Rician fading with K-factor equal to K and second moment equal to 1, it can be 
shown that 

E[Ul] = a- + 2p,a^ + 

which means that cti = • 

Substituting the objective function of (13) with (14) leads to the following optimization 
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Algorithm 3: Constellation design: Perfect knowledge of the first, second, fourth moments 


[t*^V*] = Bisection( ); 


function; V = ConstellationDesign( t) 

Pi = 0; dL.i = oo; = oo; 
dR,k = i/2ts{0); 


for fc = 2 , 3, • • • , L do 

Find the smallest pk > Pk-i such that 

if k L then 

dfl.fc = y^2ts{pk); 

end if 


2 ( \/s(Pfc) +) 


= t; 


dL,k = \/2ts{pk); 

d^k = [Pk + 0-2 - dL,k,Pk + 0-2 + dn^k]; 

end for 

end ConstellationDesign 

P = [pi,P2,--- ,pl]; 
return V^V; 


problem 

maximize min {iLAdL,k)jR,k{dR,k)) 

{pk^dL,kAR,k}kG[L] kG[L] 

0 < Pi < P2 < ■ ■ ■ < Pl, dL,k > 0, dn^k > 0 (15) 

k=l 

Note that the objective of problem (13) has been substituted in (20) with an expression which is 
still non-negative and non-decreasing in d for a fixed p^ and non-increasing in p^ for a fixed d; 
i.e., all the properties and arguments that led to Algorithm 2 are still valid. Thus, the approach 
of solving this problem is similar to the one presented in Section IV-B, with the only difference 
being that both IpAdiA and lR,k{dR,k) exhibit an easily-interpretable dependence on p^ and 
dL,k,dR,k- Algorithm (3) contains the simplified algorithm and Appendix B shows the detailed 
proof. 

Note that the suggested design leads to an algorithm which can be employed in very general 
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channel models even when there is no knowledge of the underlying channel distribution, and 
that the approximation gets better as L increases. This is because, as L increases, the transmitted 
powers are packed closer together, the decoding regions (i.e., [dR^k, di^k}) get smaller, and thus 
the approximation (14) gets tighter. Figure 3-(b) shows this numerically with Vk,^ and Vk^ as 
a function of the constellation size L in Rician fading with iF = 0 dB and 7 = 5 dB. We see 
that, with increasing constellation sizes, the approximation (14) gets tighter, which means that 
both designs lead to similar error exponents, and to approximately the same SER. 

D. Robust constellation design 

Assuming perfect knowledge of the channel distribution or even its first few moments may 
not be realistic due to changing propagation environments associated with user mobility and 
estimation errors. This motivates the need for designs which take into account uncertainties in 
the channel statistics. We build upon the design principles laid out in the previous sections to 
develop a design that performs well even in the face of channel uncertainties. 

To exemplify the suggested approach, we are going to assume that the terminals have estimated 
the first few moments of the underlying channel distribution up to some bounded uncertainty. 
We want to find a constellation that could work well for all channels with moments inside their 
uncertainty region. Recall that 

E \Ul] = s{pk) = aipl + a2Pk + ^ 3 , 

where 02 = 2(T^ and 03 = cr^. Thus, for a fixed pk, E[(7|], and hence the rate function 
approximation depends on the channel and noise statistics only through ai and a. We then 
define the following set 

J~ {(ttljCr) . Omin ^ Oil ^ Clmax) ^min CT < CTmax}) 
and note that for each / = (di, d) G we can define 

Sf{p) =di/ + d2P + d3, 

where 012 = 2 a-^, 0:3 = 
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Algorithm 4: Constellation design: Robust constellation design 


[t*^V*] = Bisection( ); 


function; V = ConstellationDesign( t) 

Pi = 0; Co = —oo; cl = oo; 
for fc = 1,2, • • • , L — 1 do 

Cfe = sup^gjr + ti/sfipk)^ +Pk; 

Find the smallest pk+i > Pk such that pk+i — Ck — sup^gjr (^t*sf (pk+i) — > 0; 

^k - [Cfc—1 5 C/i;] , 

end for 

II = [cl-1, Cl]; 

'P = [Pi,P2, ■ ■ ■ ,Pl]; 
return V = ,Il]; 

end ConstellationDesign 


Then, in order to maximize the approximate worst-case rate function for all possible channels 
inside the uncertainty region, we modify problem (20) in the following way: 


maximize 

{Pk,dL,k,dR,k}kelL] 


min 

/eJ-,fce[L] 


'^L,k \ 

sAvkY Sf{j)k) ) 


0 < Pi < P2 < ■ • ■ < Pl, d^k > 0, dn^k > 0 (16) 

k=l 

In Appendix C we show how to solve problem (16) and design a constellation which maximizes 
the error exponent for all statistics in The main difference between this design as compared to 
the previous two algorithms is the need for using power levels and decoding regions which would 
work well for any channel statistics inside the bounded uncertainty of the channel’s moments. 
To satisfy this, the consecutive power levels and decoding regions are generally spread out as 
far apart as the worst channel requires. Note also that if there is no uncertainty. Algorithm 4 
reduces to Algorithm 3. An important aspect of this approach is that problem (16), in contrast 
to the problems (13) and (17), may be infeasible. Such an example is presented below. 
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1) Infeasibility of the robust constellation design problem : In this section we present a simple 
example that shows that, for a fixed average power constraint P, a very high uncertainty on the 
channel statistics could lead to infeasibility in the robust constellation design problem (Section 
IV-D). Consider the case of constructing a constellation with L = 2, an uncertainty region 
G (e, for some e > 0 and perfectly known ai = 1 for simplicity (the case of Rayleigh 

fading). Fix t* > 0. Then, based on Algorithm 4 we choose pi = 0 and ci = \ -\- ^. We next 

choose p 2 to be the smallest p > 0 that satisfies 

p——^ -sup (t*\/p‘^ + 2a^p + > 0 p — -sup f*p + p(t* — 1)) > 0. 

If t* > 1, then p > t*p + 2^ which is impossible, and if t* < 1, then p > ^ — e. Thus, the 

smallest choice of p 2 that can be chosen is p 2 = \ — e. In this case, since 

1 + 1 

Pi+P 2<‘2P ^ Y—fff - e < 2P, 

it follows that, no matter how small t* is, if the uncertainty is so large such that 2P < ^ — e, 
the robust design problem will be infeasible. 

V. Numerical examples 

This section contains simulation studies which demonstrate and compare the performance of all 
the constellation designs proposed in this work. In the following sections, for the constellation 
designs which depend on an underlying channel, i.e., 'Pir, 7 , the {K,y) channel is 

referred to as the nominal channel, whereas any other channel is referred to as a mismatched 
channel. 

A. Comparison with a pilot-based system with PAM and a noncoherent system with ASK 

Consider a block-fading Rician fading channel hi ~ CJ\f{p, a^) with coherence time T and 
with n antennas at the receiver. We assume that both the transmitter and the receiver know the 
channel statistics but not the exact channel realization. 

In the first numerical example we compare the proposed Pir,^ with the ASK constellation, 
and with a system that uses a PAM constellation (referred to as PAM system) assuming a binary 
reflected Gray Code (BRGC) [17]. In the PAM system the transmitter uses the first p slots of 
each coherence interval to transmit pilot symbols. Based on the received signals in these slots, 
the receiver derives the MMSE channel estimates {hi} at the end of the p learning slots. Using 
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Fig. 4; (a)-(c) Minimum n for a target BER = 10 (d) Error exponents of ASK and VK,'y 


these estimates it deeodes the symbols transmitted during the remaining T — Ti slots of the 
eoherenee interval. Note that, assuming a eonstellation size of L, the effective rate of such a 
system is 

in 

The noncoherent system that uses ASK i.e., amplitudes that are equally spaced apart, performs 
decoding using an energy-based ML receiver. 

Figures 4-(a) to 4-(c) plot the minimum number of antennas needed to achieve an uncoded 
BER = 10“^ for different K and 7 = 10 dB for different coherence times T. We make the 
following observations: 

• First, our noncoherent constellation design performs significantly better than the system 
with an ASK constellation. For example, in Rayleigh fading our constellation needs 
approximately half the number of antennas to achieve the same BER performance. 

• Second, even in the case of high K (4-(b),4-(c)), which is known to all the receivers, our 
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system performs better than the noncoherent PAM system (i.e., T; = 0) which exploits the 
phase of the LOS component of the channel. Note that this is not the case with our system 
which only uses the A'—factor to decode the symbols. 

Note that Figure 4-(a) does not show the performance of the PAM system since in Rayleigh 
fading the latter cannot reach the BER target for any number of antennas as the phase of the 
transmitted symbol is completely destroyed. We also observe that, for short coherence times, 
the proposed constellation design still requires a smaller number of antennas to reach the BER 
threshold compared to the PAM system with T; = 1. On the other hand, for higher coherence 
times, the PAM system achieves better performance since the gains of learning are more than 
the corresponding decrease in the effective rate. Yet, observe that for small effective rates, e.g., 
1 — 2 bits/symbol, the additional number of antennas needed by the energy-based system to 
achieve the same BER as PAM is not more than 20. This shows that even a simple energy-based 
architecture design at the receiver, which requires only envelope detectors, could be enough to 
transmit information as reliably as a typical pilot-based system, especially in channels with small 
coherence times and high EOS, without the need for significantly more antennas. 

Eigure 4-(d) plots the error exponent Ig for different values of 7, and Rician channels with 
K = {(X), 0} dB and L = 4, for two noncoherent systems that use ASK and VK,y We observe that 
for all channel conditions, our noncoherent constellation design achieves a much higher than 
the ASK constellation. Also, for high SNR, the error exponent that uses the ASK constellation 
is not increasing as fast as the system with VK,'y This is due to the fact that the power levels are 
fixed, and do not adapt to the channel conditions. This is not the case with the Vk,'^ constellation. 

B. SER performance comparison ofVx,'^, ASK 

In the second numerical example (Eigure 5) we present, as a function of n, a numerical SER 
estimate for a 3—bit constellation (L = 8) for channels with K = —00 dB, i.e., Rayleigh fading 
and 7 = {5,10} dB. We consider the following constellations: Pir, 7 > and ASK. 

As expected, Px ,7 achieves better SER performance than all the remaining designs. Yet, the 
difference of the approximate design Vk,'^ from Vk,'^ is not significant, especially at low SNR. 
Also, the minimum distance design is significantly worse than any other design, except 

for very low SNRs, where the gap in the performance is smaller. 
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Fig. 5; SER performance comparison of VK.-y, 'pR.-y, "P™™, ASK. 


C. Performance of the robust constellation designs on the nominal and mismatched channels 

In the third numerical example we demonstrate the inefficiency of the constellation in 
a mismatched channel and the ability of to sustain good performance. Specifically, we 
consider L = 8, the case of a user with 2—dB of uncertainty in both K and 7 values and 
that the center of the uncertainty interval corresponds to the (—10,10) channel. This channel 
is approximately Rayleigh fading (K is very low) with a high SNR value. In Figures 6 -(a), 
6 -(b), 6 -(c), 6 -(d) we plot the Monte Carlo SER estimate of the P- 10,10 and P_io,io designs on 
the {(—9, 9), (—9,11), (—11, 9), (—11,11)} channels respectively. Observe the huge performance 
loss that occurs due to the overestimation of the SNR. Smaller performance loss is observed due 
to the uncertainty on the value of K, or when the SNR is underestimated. 

Figure 7-(a) presents the SER performance of Pi^io.io ^’^d P_io,io in the (—10,10) channel to 
show that even with nominal statistics, the performance of the robust design is close to that of 
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Fig. 6: SER performance of the robust constellation designs in mismatched channel 


the design that is explieitly optimized for the nominal statistics. This shows that the maximum 
performance loss due to the robust design compared to a constellation optimized for a known 
channel is tolerable, especially considering the fact that not taking into account the uncertainty 
could lead to a significant performance deterioration as presented in the previous numerical 
example. 

D. Performance on a Nakagami-m fading channel 

We now show an example in which using Vk,'^ designed for a Rician fading channel leads 
to a worse performance compared to a Vr,-^ in a Nakagami-m fading channel. This shows 
that not taking into account the uncertainty in the channel distribution, and over-optimizing the 
constellations for the Rician channel, could lead to worse performance than a much simpler 
constellation design which is based only on the first four moments of the channel. Specifically, 
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consider the ease of a ehannel for whieh it holds that 

This ehannel eould eorrespond to a Rieian ehannel, i.e., hi ~ CM , or a 

Nakagami-m ehannel with = 1 and m sueh that = \J Figure 7-(d) plots 

Ps for 7 = 0 dB and N = 6 dB in a Nakagami-m fading ehannel using L = 8, for the following 
two seenarios: 

1) a Rieian fading ehannel model and the Pg.o eonstellation design 

2) Only the first four moments are perfeetly estimated and used in the design. 

Observe that assuming Rieian fading and using the eorresponding eonstellation, leads to a worse 
performanee in Nakagami-m fading than using a eonstellation design whieh takes into aeeount 
only the first few moments of the ehannel distribution. 
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l + K 
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VI. Conclusions and Future Work 

We have formulated and solved the single-shot constellation design problem for a noncoherent 
SIMO system with a large number of antennas and an average energy-detection-based receiver. 
We present asymptotically optimal constellation designs with respect to the achieved error 
exponent when the system has perfect knowledge of the channel statistics. Then, we present 
a constellation design which requires only the knowledge of the first four moments of the fading 
statistics. The gap to optimality with this design becomes smaller with larger constellation sizes. 
Lastly, we present a robust counterpart of our designs which takes into account the uncertainty 
in the channel statistics. We exemplify the performance of all the proposed constellations, and 
compare them with existing symbol-by-symbol noncoherent schemes in typical scenarios as well 
as in pilot-based schemes. The results show that our designs are better than the noncoherent ASK 
modulation schemes and that they exhibit better BER performance than a pilot-based scheme 
that uses PAM modulation when the coherence time is small. The proposed system asks for a 
very simple encoding and decoding and for a receiver which only measures the average received 
energy across all antennas. 

Our findings here suggest that simple receiver architectures are promising alternatives to 
complex coherent designs for the large antenna systems of the not-too-distant future. We did not 
however explore the full range of optimizations that could potentially be carried out in such a 
setup. We list some directions for future research in the following: 

• Antenna correlation and how this affects the performance. This is especially relevant as 
antenna form factors go down with increasing numbers of antennas. 

• Constellation designs for a multiuser noncoherent SIMO system. Initial results towards this 
direction appear in [21]. 


Appendix A 

To begin, without loss of generality, index the constellation points such that 0<pi < P2 < ••• < Pl- Then, fix 
any codebook V which satisfies the average power constraint, and solve (13) over only the decoding regions {2i}, 
i.e., over {dL,k,dR,k}k=i- This subproblem can be written as 

maximize min {lL,k{dL,k),lR,k{dR,k)) 

{dL,k^dR,k}kG[L] k£[L] 

subject to dL,k+i + dn^k = Pk+i - Pk, Vfc € [L - 1], 

dL,k > 0,dR^k >0, Vfc e [L], 
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Observe that (17) is separable to L — 1 optimization problems, one for every k = [L — 1], since for each k, 
the constraints are separable. To see this, the constraint for k is dL,k+i + = Pk+i — Pk and for fc + 1 is 

dL,k+2 + dii^k+1 = Pk+2 — Pk+i', the former is a linear constraint between dL.fc+i, and the latter another 
constraint between dk,k+2, dfi^k+i- Each one of the resulting subproblems identifies the boundary between the 
and pk+i constellation point such that the minimum between those two is maximum^: 

maximize min{lR^k{dR,k), lL,k+i{dL,k+i)) 

dL,k + l, dji^k 

subject to dL,k+i + dn^k = Pk+i - Pk, 

dk.k+i > 0 , dfi^k > 0 . 

Each of the above problems is solved for dk^k+i, dn^k such that 

lR,k{dR,k) = lL,k+l{dL,k+l) ^ lR,k{dR,k) = lL,k+liPk+l — Pk — dR^k)- (19) 


Note that such dR^k always exists since, for any pk < Pk+i, lR,k{d) and lL,k+iiPk+i — Pk — d) are increasing and 
decreasing functions of d, respectively, with = 0, II,fc+i(0) = 0 and lL^k+i{Pk+i - Pk) > 0, lR,k{Pk+i - 

Pk) > 0. In other words, increasing dR^k increases the right error exponent of the power level, but decreases 
the left error exponent of the (fc + 1)*^ level, and there always exists a dR^k for which both are equal. Therefore, 
for any fixed V, the best decoding regions between any consecutive constellation points can be calculated by (19). 

Then, the following optimization problem finds the optimal V: 


maximize t 

{Pk}kG[L] ,{dR,k}k€[L-l] ,t 


subject to lR,k{dR,k) = lL,k+l{Pk+l — Pk — dR^k), Vfc e [L — 1], 

lR,k(.dR^k) >t, V/c G [L — 1], 

1 ^ 

^ ^ 0 < Pfc < Pk+1, 0 < dn^k < Pk+1 -Pk, yk g[L- 1]. 

^ k^l 

The solution of (20) corresponds the largest t* such that the following problem is feasible: 
find {Pk}k=i,{dR,k}kZi 


( 20 ) 


subject to lR,k{dR,k) = lL,k+iiPk+i — Pk — dR^k), Vfc G [L — 1], 
lRAdR,k)>t*,\/kG[L-l], 

1 ^ 

^ ^Pk = 1, 0<Pk < PkAi^ 0 < ^R^k < PkAi -Pk, VA; G [L - 1]. 

^ k^l 

Observe that for f* = 0 the above problem is always feasible since Ir, kid) > 0 and Ik^kid) > 0. Also observe 
that for f* = oo it is infeasible due to the finite power constraint and the fact that iR^kid), lL,kid) are increasing 
functions of d. 


'assuming dL,i = oo and dn^L = oo. 
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The problem now is to find the largest t* for which (21) is feasible. We are going to describe in detail the 
algorithm that finds whether problem (21) has a feasible solution for any fixed and finite t* > 0. The basic idea 
of this construction is that, for any fixed t*, we should find the constellation with the smallest average power 
constraint, as this is the only constraint that could lead to infeasibility of (21). To start, fix > 0 and choose 
p\ = 0. Choosing a higher value for pi can only make the problem more difficult since i((i) is decreasing in pi, 
and p 2 > Pi, thus also the rate functions for p 2 and the rest constellation points will be lower. Then, find cIr^i > 0, 
denoted as j, such that 

lRAdR,i)=t*. (22) 

The above equation has always only one solution since is an increasing function of d, /_r,i( 0) = 0 and 

t* > 0. The solution of (22) leads to the closest point to the right of p* = 0 which should be used as a boundary 
point for the first constellation point. Using a smaller boundary point would lead to a smaller rate function than t* 
since lR i{d) is increasing on d. Until now we have specified p*, j. Now, we find the smallest P 2 > P* + d*j^ 
denoted as P 2 , such that 

lLAP2-p*i-d*R^,)=t*. (23) 

Note that for P 2 = P* + i, .?l, 2(0) = 0. If there is no p 2 > P* + d*^ that solves (23), then problem (21) is 

infeasible for this t* and we need to repeat the construction for a larger t*. Note that if this is the case, i.e., if 
Il, 2 {P 2 ~ Pi~ d*Ri) < t* for all p 2 > Pi + then choosing a d’^> d*j^ ^ in the previous step would not have 
made (23) feasible. This is because, for any fixed P 2 > P* + d*j^ i, 

t* > IlAp2 -pI- rffl.i) > IlAp2 -pI- 4.i): (24) 

since lL.2{d) is an increasing function of d. On the other hand, if (23) has a solution, we use P 2 and d*j^ i to find 
d*L ,2 hy d*j^ 2 = P2 ~ Pi ~ d*R i- 

We calculate the remaining {p ^}^^3 and decoding regions iteratively by first finding djj > 0 that solves 
lR,k{dR,k) = t*, and then finding the smallest pl_^_-^ > Pk + dlj^ which satisfies lL,k+iiPk+i ~ Pk ~ d*R^k) = 

By construction, this solution corresponds to the constellation points with the minimum sum power which achieves 
a minimum left and right rate functions at each constellation point of at least t*. 

In words, this is the case since choosing a smaller dR^k than c?^ or a smaller pfe+i than p%_^_i, would lead to 
a smaller right error exponent than t* for the point, or a smaller left error exponent than t* for the (fc + 1)*^ 
point. Then, if it holds that ^ J2k=iPk — problem (21) is feasible for t*. Identifying the largest t* for which 
(21) is feasible solves (13). 

To efficiently perform this procedure we can employ a simple bisection algorithm (Algorithm 2). To see this, 
observe that for a i, such that t* < i, the corresponding constellation design leads to higher (or equal) average 
transmitted power (infinite power if the problem is infeasible). This is true because, at each step of the constellation 
design, finding dR^k that satisfies lR^k{dR^k) = t will lead to a with > d*j^^, and finding pk which 
satisfies lL,k{Pk - Pk-i - dR,k-i) = i will lead to a pk with Pk > pI- □ 
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Appendix B 

In this appendix we take into account the approximation of the left and right rate functions shown in (14) 
to simplify the algorithm needed for a constellation design that uses only the first, second and fourth moments. 
Equation (19) can now be written as follows 

_ _ {Pk+1 — PkY _ 

which means that the feasibility problem (21) is simplified to 

find {Pk)k=i 

subject to pfe+i -Pk> (^^/s{pk+i) + \Js{pk)^ ^ 25 ) 

k=l 

Then the procedure described in Appendix A is now simplified to the following: Fix t* and choose p* = 0. Then, 
iteratively choose the smallest pfc+i > p^. for A: = 1, 2 • • • , L — 1, such that 

jPk+i -Plf ^ 

2 {ys{pk+i) + v^iPk)) 

If no pk+i > P*k exists, then problem (25) is infeasible. □ 


Appendix C 


In this appendix, we show the details of the robust constellation design problem. This problem is simplified if we 
denote a constellation using {pk}k^i and {ck}kZi^ where Cfc is the boundary of the decoding region between the 
Pk and the pfc+i constellation point. Then, using the approximation shown in (14), the problem of maximizing the 
worst case approximate rate functions for all the channels inside the uncertainty region T is expressed as follows: 


maximize t 


subject to 


Ck-Pk> sup + t^Sfipk)^ 
Pk+i - Ck> sup tJsf{pk+i) - 

V ^ 


,Vfce [L], 

aA ,Mk e [i- 1], 


1 ^ 

j^'^Pk = i.p/c > o,Vfc e [L]. 


(26) 
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(27) 


This problem is equivalent to finding the largest t* > 0 which gives a feasible point in this formulation; 
find {Pk\i^i,{ck}lzl 

subject to Ck — PkZ sup ( + t*\/sf{pk) ) ,VA: € [L], 

^ ) 

Pfc+i - Cfc > sup ^ Sf{pk+i) - ,VA: G [L- 1], 

1 ^ 

J^'^Pk = ^,Pk > 0,Vfc G [L]. 

fc=l 

Solving the above feasibility problem can be done as follows: Fix a small f* > 0 and choose pi = 0 and d*j^ ^ = oo 
so that Ii{d*j^ i) = oo. Using p* > 0 would lead to a sub-optimal solution since the transmitter has an average 
power constraint and s / (p) is an increasing function of p for every / G -F. Then, choose c* which satisfies 


4 = sup + Sf{p\)^ +pI, 


and as pj, the minimum p that satisfies 


p - c* - sup \Jsf(p) - > 0. 


(28) 


Note that for 0 < t* < inf --d= there always exists a p > — infthat satisfies the above equation. To 

see this, define the following auxiliary function Wf(p) = p — c\—t* \/^fip) + , for which, for any fixed / G 

it holds that Wf{cl — a^) < 0 and lim = x _ > 0. 

Note that choosing a higher value for c* would only make p^ larger (or infinity) and thus, use more transmit 
power than necessary (or make the problem infeasible). Using the same procedure we can sequentially specify all 
{PfciLi and Then, if 

k 

the problem is feasible. However, if the average power constraint is not satisfied, it is not possible to guarantee 
this error exponent for all channels in J^, since in our construction, we pack the decoding regions and constellation 
points as closely as possible. To see this, if in the above construction we choose any value Ck > 4, then the 
corresponding p which satisfies (28) would be larger than pl_^_i since 

p-Ck- sup {t*\Jsf{p) - <p-c*k- sup {t*\J~s4j})- ,Vp > Pfc + Ck- 

□ 
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